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O ' Abstract. In the context of stochastic volatility models, we study representation formulas in 
terms of expectations for the power series' coefficients associated to the call price-function. As 
in Antonelli and Scarlatti [5] the expansion is done w.r.t. the correlation between the noises 
driving the underlying asset price process and the volatility process. We first obtain expressions 
fx. ' for the power series' coefficients from the generalized Hull and White formula obtained in Alos 
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[2]. Afterwards, we provide representations turning out from the approach for the sensitivity 
problem tackled by Malliavin calculus techniques, and these allow to handle not only vanilla 
options. Finally, we show the numerical performance of the associated Monte Carlo estimators 
for several stochastic volatility models. 
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§ Introduction 

^ I In the last three decades a lot of works have been produced in order to generalize the seminal 
model by Black and Scholes [IT]. The purpose of many researchers was to formalize models 
with a more general stock price dynamics able to explain the features observed in the financial 
market. Indeed, with the assumption of constant volatility (as in the Black-Scholes model) one 
^ ! cannot explain well known stylized facts like the so called smile effect, volatility clustering, fat 
tails or excess kurtosis in the log-returns' distribution. It follows that the empirical observation 
of such phenomena leads to model the volatility itself as a stochastic process with a dynamics 
correlated with that of the risky asset. The most popular stochastic volatility models are, 
among others, the models by Hull and White [IT], Stein and Stein [19] or Heston [16]. We also 
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refer e.g. to Frey [H], Ghysels, Harvey and Renault [15] or Fouque, Papanicolaou and Sircar 
[T2~| for general surveys. 

In the context of stochastic volatility models, several recent papers deal with the behavior of 
the price of derivatives w.r.t. a small parameter. For example, Benhamou, Gobet and Miri [TO] 
derive analytical formulas for pricing vanilla options for a time dependent Heston model by 
using a small volatility of volatility expansion. Another parameter of interest is the correlation 
between the noises driving the underlying asset price process and the volatility process. From 
the analysis of financial data it is in fact evident that there exists a negative instantaneous 
correlation between the volatility and the stock price process. This is exactly what we are 
interested in: we study here power series representations for European option prices. 
We assume (see Section [2] for mathematical details) that the risky asset S and the volatility 
process v fulfill the dynamics 

p(v e )d6 + J r]{v B )dB l e 
S s = e x + J rS e dO + j f(v e )S e (p dB] + v 7 ! - P 2 dB 2 e ) 

where B is a 2-dimensional Brownian motion and p G (—1, 1) denotes the correlation between 
the Brownian motion B 1 and the Brownian motion Z = pB 1 + \J\ — p 2 B 2 . In this setting, the 
well known formula for the call price function is given by 

u(t,x,v;p) = E (e~ r ^ (S^' v - K)\ . 

Antonelli and Scarlatti proved in [5] that under suitable regularity assumptions on /, r\ and p 
(that we will recall in Section [3]), u is C°° as a function of p in a neighborhood of p = 0, where 
it can be developed in power series: there exists R e (0, 1) such that 

u(t,x,v\p) = ^2g k (t,x,v) p k , \p\ < R. 

k>0 

The power series coefficients 

g k (t,x,v) = yd k p u(t,x,v;p)\ p=0 

solve suitable parabolic partial differential equations and thanks to Feynman-Kac type formulas, 
they can be represented as expectations of functionals of the diffusion pair (S, v) evaluated at 
p = 0. 

Starting from this fact, in this paper we study alternative probabilistic representations, as an 
expectation, for the Taylor coefficients ^'s. 

We first deal with the results in Alos [2], that concern a generalization of the classical Hull and 
White formula [IT] for European option prices. Starting from the Alos formula, we first give 
an alternative representation for call price function u and then we perform derivatives, getting 
an expression for the Taylor coefficients with a direct approach which leads to simple formulas. 
But from the financial point of view, it is also clear that the Taylor coefficients g^'s can be 
thought as sensitivities of the price-function u(t,x,v;p) with respect to the correlation coeffi- 
cient p. Thus, we obtain probabilistic representations for the g^s in terms of expectations by 
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means of Malliavin calculus techniques. And this brings to further alternative representation 
formulas. Here, we get the results by using the Malliavin calculus in the direction of the noise 
B 2 . This allows one to find weights which are easy to handle and worth to set up plain Monte 
Carlo estimators for the power series coefficients. It is worth to be said that one could actually 
find results from the use of both the noises B 1 and B 2 but this procedure brings to difficult 
formulas, unfeasible to be used in practice. Moreover, the Malliavin approach does not depend 
on the particular form of the payoff function, although it must depend on the asset price at 
maturity, and this is a basic difference with the result in Antonelli and Scarlatti [5] (even if 
their result can be generalized to other situations, see Antonelli, Ramponi and Scarlatti [6]). 
So, we numerically study the behavior of the Monte Carlo estimators to price call options by 
using the Taylor expansion up to the first and the second order, by means of the representation 
formulas for the coefficients we provide. It is worth to be said that all representations work 
efficiently, the Malliavin weights being the ones giving results a little bit less accurate but, at 
the same time, providing the most flexible approach. 

Plan of the paper. Section [2] is devoted to set up the model, the notations and the main 
requirements we need and Section [3] contains the results obtained by Antonelli and Scarlatti 
[5]. In Section |4] we study the series expansion from the Alos [2] generalized Hull and White 
formula whereas the aim of Section [5] is to obtain an expression for the coefficient of the 
series expansion by using the Malliavin calculus approach for the sensitivity problem. Finally, 
numerical results for several well known stochastic volatility models are discussed in Section El 

2 The framework 

We resume here the notations, the model and the requirements we need throughout this paper. 

Let T > be a fixed time horizon and let (Q, J 7 , IP) denote a complete probability space on 
which a 2-dimensional Brownian motion B is defined over [0, T], and we let (J r t)t e [o ) T] stand for 
the natural Brownian filtration augmented with P-null sets. 

We consider a market model with one risky asset and one risk- free asset (the money market), 
the latter being assumed to be modeled by means of the instantaneous interest rate r. We 
assume that the price S of the risky asset and the stochastic volatility v are modeled as follows: 
for t G [0, T], x G R and v > 0, the pair (v, S) evolves as 



for s € [t,T], where p G (—1, 1) denotes the correlation between the Brownian motion B\ and 
the Brownian motion Z t = pB\ + a/1 — p 2 B 2 . Thus, ((T]) and ([2]) both say that the underlying 
probability measure P is assumed to be a risk-neutral one. 
Let us consider the following assumption: 

(H) /, p,rj : R — >■ R are C°° functions whose derivatives of order > 1 are bounded. 




(1) 



(2) 
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In particular, (H) ensures the existence and uniqueness of a strong solution of ([I])-© and the 

processes £ s = In S s solves 



f (r ~ \f\v u ))du + J S f(v u )(pdBl + y/T=fdBl). (3) 

In what follows, it will be useful to give a special denomination to the process £ whenever p = 0. 
Therefore, we set 

is = Cs\ P =o : is = x + J" (r - l -f(v u ))du + J" f(v u )dB 2 u . (4) 

We will use the notation £l' x ' v , Cl' x ' v an d v l' v to stress the dependence of £ s , £ s and u a w.r.t. 
the starting instant and position, and this will be done for all the diffusion processes we are 
going to define. For further use, let us also recall that under (H) the first variation process 
Y s = d v vl ,v and its inverse Z s = Y" 1 , s >t, are both well defined and solve the equations 

Y ' = 1 + f **W0*W« + f v'(rt v )YudBl (5) 
Z s = l+ [ S [-»' + (rjTM v )Z u du- f \'{v^)Z u dBl 



respectively. Equivalently, 

Y s = exp ( ^ [p' - l^) 2 ] (v^)du + ^ T,'(v?)dBl) . 

We finally recall that, under (H), all above processes £ > t < x < v ; v t,v , Y, Z can be collected in a mul- 
tidimensional diffusion process X t,x,v , whose diffusion coefficient and drift are either Lipschitz 
continuous and with sublinear growth. Thus, in particular, standard LP estimates hold, as the 
Burkholder one: for every p > 1, t G [0, T] and i,f)6l there exists a positive constant C such 
that 

E( sup \X l s ' x ' v \ p ) <C{T- t) p/2 {l + \x\ p + \v\ p ). 

t<s<T 

3 Antonelli and Scarlatti power series expansion 

We recall here the result obtained by Antonelli and Scarlatti in [S] . They proved a development 
in power series of the European call option price-function with respect to the correlation between 
the noise of the risky asset and the noise associated to the stochastic volatility. 
Consider a European call option written on our market model. For a fixed maturity T and 
strike price K, its price-function is given by 

u(t, x, v;p)=E (e- r{ ?^ (e^" v - k) \ (6) 



When the underlying asset price process evolves following the Black and Scholes [TT] model, 
06]) has a closed form solution, following the worldwide famous formula 

BS(t, x, a) = e x N{d 1 {t, x, a)) - A'e"^^ N(d 2 (t, x, a)), (7) 
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where 

, x-\nK+ (r - a 2 /2)(T -t) , , , = 

d 2 {t,x,a) = \ 1 ,y - '- and d 1 = d 2 + a^/T - t. (8) 

o\/T — t 

In order to tackle the problem of the regularity of the price-function u given in (EJ w.r.t. p, the 
following assumption is considered in [5]: 

(AS) For every t G [0, T) and u£l i/iere exists a constant C± > snc/i i/iai /or any g > 1, 

v((J t T \fM' v )\ 2 dry q ) <Cf. 

Then, the following result is stated in Theorem 2.2 of Antonelli and Scarlatti [5]: 

Theorem 3.1. Assume that (H) and (AS) both hold. Let u = u(t,x,v; p) denote the price- 
function (0|). Then, 

(i) the map p t— > u(t,x,v; p) is infinitely differentiable in a neighborhood of p = 0; 

(ii) there exists R G (0, 1) such that for \p\ < R one has 



u(t,x,v;p) = ^2 —dpU(t,x,v;p)\ p=0 p 



k>0 

For k > 0, let 

g k (t,x,v) := — d k p u{t,x,v;p)\ p=0 

denote the coefficients of the Taylor expansion of u. In [5j, it is shown that the g^'s solve suitable 
parabolic partial differential equations and, by using Feynman-Kac type formulas, each g k can 
be represented in terms of an expectation of a suitable functional of the diffusion process (£, v) 
(see (3.1) in [5]). Moreover, by conditioning w.r.t. the noise B 1 , such expectations reduce to 
ones of functionals of the volatility process v only. Here, we write down the formulas for g k as 
k = 0,1,2 (for the general formulation, see (3.3) in [5]): one has 



g (t,x,v) =E | BS(t,x,\l ^^ ) | (9) 



9i 





t,v \ 
C [i,T] \ 




(M t ' v )x / 



g 2 (t,x,tO ^e-P--') / T E( 3 fatff)" 1 * (11) 

^ V((M^) [t , r] )^ V 7 



x 



da? 



jy( 5 )(d 2 ) _ iv( 4 )(d 2 ) \ / ^ /(M^) [t , 



(^' u )[t,T] VWm/V V T "* 
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where N denotes the cumulative distribution function associated with the standard Gaussian 
law and N^(-) is its nth derivative, d 2 is given in (jSJ) and for s E [t,T], 



{M t,v ) 



t,v 

\s,T] 

[«.n 



(JV)W) / U f){vf)d v vf (d^y'dpda, 

J a 



(12) 
(13) 



4 Series expansion A - from the Alos formula 

We deal here with the price-formula obtained by Alos in [2] and we derivate it in order to state 
representation formulas for the series expansion. 

Consider a European option with maturity T > and payoff /i(£r)> whose price, as seen at 
time t < T, is given by 

^ = 1(6-^-^(^)1^). (14) 

We denote through BSh{t, x\ a) the price- function of this option under the Black and Scholes 
model with constant volatility a, that is 

BS h (t,x;a) = E^^h^^-^-^^-™^ , (15) 

W denoting a standard Brownian motion on R. 

In order to state the starting result from Alos [2], one needs some requirements (see (H1)-(H4) 
therein, at page 357). In our context, which takes into account Assumption (H), we can rewrite 
them as follows: 

(P) the payoff function h : R — > R + is continuous and Cl-piecewise; 
(A) there exists a > such that \f\ 2 > a. 

Let us remark that in Alos [2] the model is more general (the volatility is not necessary a dif- 
fusion process, neither Markovian, but definitively measurable w.r.t. the noise B l ). Moreover, 
a further request appears (see hypotheses (H4) therein, at page 357), which is not written in 
above (P) and (A) because it can be actually dropped (see e.g. the more general result - jumps 
are allowed - in Alos, Leon and Vives [1], Theorem 4.2 therein). 

Theorem 3 in Alos [2] states the following generalization of the Hull and White formula [TT] . 

Theorem 4.1. Assume that (H), (P) and (A) hold. Let V and BSh be given by [Lfy and 
(T73J) respectively. Then, for any t e [0, T], one has 
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where, as s G [t, T], 

(M)[s,t\= [ {f{vr)) 2 dr (17) 

H(s, x- a) = (81 - d 2 x )BS h (s, x; a) (18) 

K=( T D]{{f{v r )Y)dr\f{v s ) (19) 



The notation D 1 in ( TT9|) stands for the Malliavin derivative in the direction of B 1 . We do not 
enter here in the details of Malliavin calculus (we will do this in Section [5]). 
Let us now specify to our problem. First, we notice that (A) implies that (AS) holds. And 
obviously, (P) is fulfilled by the call option payoff. Then, we can proceed to study the Taylor 
series expansion of the call price-function u by developing formula (1161) w.r.t. p. To this 
purpose, let us recall the notation BS for the call option price-function given by the Black and 
Scholes formula, see (E])- We first state the following result. 

Proposition 4.2. Assume that (H) and (A) hold. Then the European call option price- 
function u given in (0|) can be written as 

u(t, x, v- p) = E(rf' v ) + P n^T'\p)) (20) 

where 1 { ana 1 2 {p) are defined as 



rr> v = BS{t,x,d {M !' V)[t > T] ') (21) 



T-t 

"T j.t,v 



Y^{p) = -Ke-^ I *'\ [D^N'(D^)} (s; p)ds 



(M^) [t)T ]-p 2 (M^) [t „ 
{M t,v )[ Si T], s G [t,T], being given in / TJ^) and for s G [t,T], 



(22) 



tf v = (fvM' v ) [ (ff'MW/ (d v vYY l dr, (23) 



(M t>v ) \ 2 T] d 2 (t, x, V^P) + p J M")dBl 

D^is; p) = [ -—^ V - i m (24) 

((Af^) [t , T] -p 2 (MM) M ) 1/2 

The proof of Proposition 14.21 is an immediate consequence of the following result. 
Lemma 4.3. Under assumptions (H) and (A) the following statements hold: 
(i) the European call option price-function u given in |Qj can be written as 

u(t,x,v;p) = E(Lf ' u ) + P E(f t 2 ^'(p)) (25) 
where y\ x,v is given in l[21\) and 

(M''"}^^ and ipl ,v being given in (T73|) and respectively; 



ii) if T 2 ,x,v and T 2 ,x,v are as in j[26)) and [2^1 respectively, one has 



E(r^(p))=E(rr' v (p)) 



t,X,V i 



(27) 



Proof, (i) As already observed, requirements (H) and (A) imply that both Theorem 13.11 and 
Theorem 14.11 hold. So we can actually proceed with the proof. For a European call option, 
the price-function u is given through u(t, £ t , v t ; p) = V t , where V t is as in formula (fTB"j) with 
h(£) = (e ? - K) + . We firstly notice that from (TIT]) and CE} one has {M) [s>T] = (M*^)^. So, 
by using the Markov property one immediately gets 



E BSi t,e 



(M) 



T-t 

Therefore, we prove that for any s G (t, T] 

I A 



E BS Lx 



(Af*>«) 



[t,T] 



T-t 



x = it 

V = v t 



-2E Ke- r ( T -V 



(M) [S ,T] 

T-s 
1 



[d 2 N'(d 2 )}[s,e s ' x ' v - 



T-s 



1> 



s 



(28) 



x = ft 
V = v t 



where H and A are given by f fTSl (with BSh replaced by BS) and ffHIj) respectively. And 
moreover, we prove also that one can actually use the Fubini theorem, so that (1261) will hold. 
By recalling ((?]) and by using the well known properties d x BS(t,x,a) = e x N(di), N"(r) = 
—tN'(t) and N'(di) = N' (d 2 )Ke~ x ~ r ^ T ~ t \ straightforward computations give 

e -r(T-t) 

(8g - BS-(t , x, a) = -K 2lrr ■ [d 2 N\d 2 )} (t, x, a) 



so that 



(M) [S ,T] 

T-s 



-K 



a 2 (T-t) 

e -r{T-s) 



(M)[s,T] 



■[d 2 N'(d 2 )][s,i 



(M)[.,T\ 

T-s 



Let us consider A as in f JT9|) . First, thanks to (H), for r > u>t one has 



D 1 J 2 (v r ) = 2(ff)(v r )D 1 s v r and D]v r = (d^f (d v v^Y 
d v vl ,v , s >t, denoting the first variation process of v t,v , see fl5]). Therefore 



As = (fv)(v s ) / 2(ff')(v r )(d v vl' v {d v vT) 1 



dr = 2-0 



t,Vt 



v=v t 



the process ijj t,v being given by (12"3"|) . By inserting the expressions for i7(s, £ s ; , 
we get 



W[s,T] 



and A, 



(M) [S ,T] 



T — s 



A. 



-2AT 



-r(T-t) 



<M^) [s , r] 



■[d 2 iV'(d 2 )] s,£ 



T-s 



and the Markov property allows to conclude that (1281) really holds. It remains to prove that 
the Fubini theorem can be applied. So, we show now that for every p > 1, the above process 
Z belong to L p w.r.t. the product measure dP x ds. 

Since \xN'(x)\ = \N^ 2 '(x)\ and sup^ \xN'(x)\ < const (from now on, const denotes a suitable 
constant, possibly varying) and since (A) gives (M t,v )[ s ^] > a(T — s), one has 



[ \Z S \ P < const ■ f — 1 - \^f\ p ds. 
Jt Jt (T-s)p Wa 1 



Then, it remains to prove that E(J t T \il) t s ,v \ p ds) < oo. Notice that in particular this gives 

that ip*' v G L p w.r.t. dP x ds. One has 

|$T < const ■ (1 + K'T) \d v v t s ' v \~ p (T - sY' 1 [ (1 + \v t ; v \ p )\d v v t r ' v \ p dr 

J s 

so that 

T -> /-T f^/, , \ t,v\p\\a v t,v\p J r 

W\ p ds< / js [ 7J Jl 7 r 1 (i + i^t^X'T^- 



By using the Burkholder inequality, we know that sup f<s<T (l + |?4 ,u | p )|c>uU*'' y |~ p G L q for all g. 
So, by using the Holder inequality, it remains to prove that for a > 1, 

t ^i^e(( j\i + i^ni^r^) a ) 1/Q < oo. 

By the Jensen inequality one gets 

r-T 



E 



(1 + |^'T)|^*Tdr) a ) ° <(T - s) 1 -^ y (1 + \v t r ' v \ p )\ a d v v t r ,v \ pa dr) 

<(T - s)e( sup (1 + KT) a \d v v* r ' v \ 

V t<r<T 



1/a 

\pa 1 



i<r<T 

<const X (T — s) 

and this gives the result. Statement (i) is now completely proved. 
(ii) By applying the Fubini theorem, we have 

where J-^ = o"(£>^ : s < T). We notice that 
d 2 (s,£ 



(29) 



SS 5 



(M^) [S , T] 



T-s 



s,T] 
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Hence, conditionally on J 7 ^ 1 , d,2 (s, \j ^ M r ^j/ ^) is a Gaussian r.v. with variance and mean 
given by 

^( S ;p) 2 = (l-p 2 )- (MM)M 



(M"*'%,T] 



respectively Now, if Z ~ A/"(/i, a 2 ) , straightforward computations give 



e 2(i+ CT ^y 



(l + a 2 )i 

Thus, by inserting a = o~ t,v (s; p) and p = p t,v (s] p), one gets 



E([d 2 Ar'(d 2 )]( S) ^-'/^^ 



1 



T-s 



T 



s J 

(M ' V )[ S ,T] 



(30) 



(M^} m -p*(M^} [t , 



where D t,x ' v (s;p) is defined in (1241) . By putting (1301) in the use of the Fubini theorem 

allows to get the claimed result. So, in order to conclude we must only to show that the Fubini 
theorem can be actually used, and this holds if 



E 



But this is true because ( M t,^ T] _p2{ M M) [t ' [-P^'^^-P*' 3 ^')]^; P) is bounded and e L p 
w.r.t. d¥ x ds. □ 

Now, under (H) and (A), Theorem 13. II can be applied, so that in particular p H- E(r 2 ' a: ' i; (p)) is 
actually C°° in a neighborhood of p = (we stress that Y t { x,v is independent of p). So, in order 
to give an alternative representation of the coefficients of the Taylor expansion, we compute 
the coefficients 

u k (t,x,v) = d k p u(t,x,v;p)\ = ^(pE(rf< u (p))Y , k > 1 (31) 

by taking the derivatives inside the expectations, and this can be done. Of course, u k = k\gk, 
the gkS as in Theorem I3.1[ althought their representations in terms of expectations are not 
generally the same. In fact, we have 

Proposition 4.4. Assume that (H) and (A) hold. Let T^'^p) be defined in (TJ/SP. Then for 
every k > one has 

W«M) |( _ = Ke-^-'>E( J ^- ] £ S?*(s) ds 
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with 



nt,X 



j=0 



in which {G\ x ] f ) }k and {G t { L k ,v }k are given by 



t.x.v^ 



(32) 



and 



£ = 0: G^(s;p) = 

e>l: G^{s,p) = Y,li[N {2+ "\d 2 )](t 



T-t 



!/=l 



T-t 



(33) 



hi > 1 Vi t=l 
Til + ... + /i„ =£ 



respectively, where 



d 2 [t,x, 



(2n)\ f (M t < v ) [ tt8 ] \ " 
(n!) 2 4« V(M^) [t!T] / V T - 

(2w)! / (M'^.n n f/^dBi 



if h = 2n, 
ifh = 2n + l. 



(34) 



As a consequence, one has 

u n (t,x,v) := d?u(t,x,v;p) lp=0 = nKe^^E^ ^ - jf 5ffi(a) <fe) , n > 1. 



Proof. One has 



/<'p-'-( t '-*) /- T 

G{ (s; p) G 2 ' (s; p) W (is 



with 



Gf{s-p) 



and Gf'^s- p) = N"(D^ v )(s; p) 



(M*>% iT] - f?{M*,% 
(recall that —tN'(t) = N"(t)). We write formally 

h 



3=0 



k 



E 



1 



j J \{M^) [t , T] J t 
So, in order to prove our result it is sufficient that if p is close to then 



d k -^G t f' v {s-p)& p G t t v {s-p)^ds 



E 



{[ \d k r J G^ v ^p)^ p G^{s-p)^\ds) < oo 
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(recall that (A) gives 1/ {M t,v )uT} < const). In the proof of Lemma [4731 we have already proved 
that E(/ \^ v \ p ds) < oo for all p. Therefore, it remains to show that there exists R > such 



that if \p\ < R then 

e( ^ ^G'f'^s] p)\ P ds} < oo and e( ^ \d 3 p G? ,v (s; p) \ P ds) < oo 



for all p. So, we prove these facts. And to simplify notations, we drop the dependence on t, x, v. 
From (1 — t)^ 1 = ^ n>0 i n , < •"•> one immediately gets 

p) = E Gi,k( s ) P k where {G lik (s)} k is given in 



fc>0 



Since sup fc sup se [ tT ] \Gi t k(s)\ < 1, we can write d p Gi(s; p) = J2 k>0 G lik (s) d p p k . So, for any 
< R < 1 and |p[ < R, 



C^r being a suitable positive constant, and this gives the first result. Let us now study G 2 (s; p). 
We can write 

l Af(2+v)( J~)( q . n \\ v 1 

d> p G 2 (s;p) = di p {N"(D)}(s;p)=j\Y, II XT^*^ ^ 

v=\ hi > 1 Vi i=1 *' 

/ll + ... + = j 

Now, for > one has N^ 1+k '(r) = (—l) k H k (r)N'(T) where denotes the Hermite polyno- 
mial of order k. So, for some positive constant Lj, it holds that sup^. |A r ^ 2+1/ - ) (a;)| < Lj for every 
v — 1, . . . , j. Therefore, 

iW«;p)i<¥Ei E n^l*^^)!- (36) 

f=l ' hi > 1 Vi i=l l ' 
hi + ... + h v = j 

Consider now D(s; p). Since for \t\ < 1 one has (1 — t 2 )~ l l 2 = J2 n >o Cn ^ n ^ w hh c n = /j)" |n ; by 
inserting in formula (I24|) for |p| < 1 one gets 



D(s; p) = E jDfc ( s ) P fc wnere {-Dfc(s)}fc is given in 

fc>0 

And since for every < R < 1 there exists C > such that |c n | < CR~ n , we obtain 



iDfcfs)! < CM s R~ k , with M, 



d 2 (t,x, 







i: /(««) ^ 




(M)[t,T) 
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Moreover, for \p\ < 1 one has d e p D(s;p) = ~^2 k>0 Dk( s ) d e p p k . So, if we take \p\ < R/2 we can 
write 

u 

\S^D{a;p)\ <Y, CM ^ R ~ k Tk^KM k ~ l ^ ^ +X CM s BT l . 

k>£ \ )■ 

As a consequence, we get 

\dlG 2 {s-p)\ <j\Lj J^^ Yl f[CM 8 2 hi+1 R-^ ■ < Cj, R (l + 2CMJ 

v=\ ' hi>lVi *=1 
hi + ... + h v = j 

where Cj t R denotes a suitable positive constant depending on j and R. Since E(/ t \M s \ p ds) 
< oo for all p, one gets the desired integrability property for G^s; p) as well. 
So, we can actually write 

9jE(r^(p)) u=o =ife-^-^( * y{j^ V - f d k ^G t r{s-p)^G t t v {s ] p)^ds) 

j — Q \ / [^>^~1 t \p — 

Now, 

d k -iG 1 (s;p) lp=0 = (k-jy.G^.jis) and 

"=1 hi > 1 Vi i=l 

hi + ... + h v = j 

so that the statement actually holds. 

Finally, for n = one has Uo(t, x, v) = u(t, x, v; 0) = ~E(T t { x ' v ) and for n > 1 

M n (t,x,w) = ^M n (t,x,w;p)| p=0 = ^""^(r^p))^ 
and the proof is completed. □ 

By using Proposition 14. 4[ we can now immediately deduce a representation for the coefficients 
Mo, U\, U2- As already observed in [5], from the numerical point of view these are the relevant 
objects which the attention has to be focalized to. 

Proposition 4.5. Suppose that (H) and (A) hold and let u k be as in l[3l\) . Then the repre- 
sentation in terms of expectation for uq and u\ both agree with the ones for go and g\, given by 
(TJ|) and IfTu)) respectively, whereas for k = 2 one has 

u 2 (t, X ,v) = iKe^E^^N^id,)^ V®^) (37) 

where 

£$ T] = J Iff £ f{vf)dB] da 
the process ip t ' v being defined in 
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Proof. Uq is trivially obtained putting p = in (|20|) . and this brings to the same representation 
of g in Regarding U\, Proposition 14.41 gives 



= G%(s) G^'(s) = [N^(d 2 )] (t, x, \\ {M ^Y ] 



and one gets u\ = g\. Moreover, 



-1 (*) = ^l,o( S ) G 2,l ( s ) = l N{ '^(t, X, \ J i /2 > 

hence formula fl^Tj) holds. □ 

We observe that the second order coefficient u 2 has a representation in terms of expectation 
which is really different from the one for g 2 , as it follows by comparing (157|) with (ITTj) . Moreover, 
we would like to stress that representation fl37|) for u 2 is drastically simpler than expression 
(fTTj) for (72: this permits us to perform numerical analysis with a much smaller effort. 



5 Series expansion M - from Malliavin calculus tech- 
niques 

This section is devoted to the representation of the derivatives in terms of expectations by using 
the sensitivity representation method due to the Malliavin calculus approach. 
We first recall some notations and results from Malliavin calculus, and we refer e.g. to Nualart 
[IB] and Bally [7] for this topic. Here, we perform the Malliavin calculus restricted to the time 
interval [t, T], that is of interest for pricing options. So, once for all, we fix the starting instant 
t G [0, T) and the starting points i,d£|. 

For i G N and p > 1, we let Tb e ' p stand for the space of the random variables which are £-times 
differentiable in the Malliavin sense in LP . For F G V> e ' p and for a multi- index a such that 
\oc\ — j < i (that is, a = (ati, . . . , <x,) G {1, 2} J - recall that we are dealing with a Brownian 
motion on M 2 ), D a F denotes the Malliavin derivative of F corresponding to the multi- index a. 
In particular, the notation DF = (D 1 F, D 2 F) stands for the first order Malliavin derivative of 
F and as i — 1, 2, D l s F, s G [t, T], gives the derivative of F in the direction of the z-th Brownian 
motion (roughly speaking, the derivative is done w.r.t. the infinitesimal increment A5*). As 
usual, we set IT/' 00 = n p >iD^ p . Sometimes, for the sake of notation, we need to define and 
D 0,o °. So, we put = L p (tt, ^ T ,P) and D 0,o ° = n p >i©°' p . 

Moreover, we let Dom 2 (S) denote the domain of the Skorohod integral, which is the adjoint 
operator of the Malliavin derivative in D 1,2 . And more generally, Dom p (5) denotes the set of 
the processes which are integrable in the Skorohod sense in LP . As usual, we put Dom^S) = 
r\ p Dom p (5). For further use, we also consider the set ll ,p of the 2-dimensional processes g G 
L 2 ([t, T] x Q) such that g s G H/' p for all s G [t, T] and for every multi-index a with \a\ = j < £ 
(that is, a G {l,2p') one has \\D a g\\ L 2 {[t!T]j +i ) G L p {Vt). Recall that L 1,p C Dom p (5). We define 

1/' 00 = r)pL e ' p . 
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For G G D^' p , ip G Dom p (§) and g G h e,p , we recall the Malliavin-Sobolev norms 

\\G\\i P =n\Gy)+Y, E E ((/ T J^ G i 2dsi '--^) P/ l' 



3=1 a : |a|=j 



p/2n 



ip=E(i^)n+E((y tr i^i a is) p ), 

\9s\ 2 dsy E ®((J tTj JK.., j 9s\ 2 dsi---ds j ds) ). 



3=1 a : |a|=j 



We also recall the following Holder type inequalities: for I G N and p > 1, there exist constants 
a^p > and > such that for every a, (3 > 1 with - + ^ = 1 one has 

ll^iG^lkp — a e,p \\Gi\\i tPa \\G2Wejpp (38) 
\\G^\\ s , p <b e>p \\G\U !pa \\<iP\\ 5iP p (39) 

for all Gi,G 2 ,G G and if; G Darned). Moreover, if G G D m >°° and g G l/>°° then 

J T (D S G, g s ) ds G D £,0 ° and for every p there exists a constant q )P such that for every a, /3 > 1 
with - + \ — 1 one has 

a p 

T 

(D s G,g s )ds < Q i p||G||/ + i iPa ||0|| L «j,/j. (40) 



Let us remark that inequality ( 140]) is perhaps less popular but it can be straightforwardly proved 
by using the same technique allowing one to get fl38|) and f )39|) . 

To our purposes, we need an integration by parts formula which is the starting point of our 
results, and it is well known that Malliavin calculus produces integration by parts formulas. 
It is not the unique one but this gives formulas for the series expansion coefficients which 
are feasible to be used in practice (see next Remark 15.61 for further details). And in order to 
proceed, we recall the notation C™(M) to denote the set of the infinitely differentiable functions 
: K — > M with compact support. But we need also a condition, slightly stronger than (AS), 
which is the following: 

(M) For every t G [0,T] and v G R there exists a constant C > such that for any q > 1, 

e(| |/(^)|-%) < c 

Notice that, under (H) and (M), one has l/f(v^ v ) G L £ '°° for all I G N (actually, the right 
writing should be g G L/'°° where g\ = and gi = 1/ f(v t,v )). Then, we have: 

Lemma 5.1. Assume that (H) and (M) hold and let t > 1. Then the following integration 
by parts formula holds: for every G G one has 

E(<p\^ v )G) = E(<p{& x > v ) H t,v [G\) V / G C °°(R) 
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where the weight T-L t,v [G} is given by 

m which <5 2 stands for the Skorohod integral over [t,T] w.r.t. the Brownian motion B 2 . More- 
over, n^ v [G] e if- 1 ' 00 . 

Proof. By the chain rule, for s £ [t, T] one has 
so that 

7 - * it f{Vs ) 

Since G £ D 1,0 ° and l/ZO^) £ .0077100(52) then G/f(y*> v ) £ Dom^^), so the duality rela- 
tionship gives 

But one has also that l//(t> < ' 1 ') £ I/' 00 and 

^(77^t)= g fjn^ dB2 s- fink*** 

which says that / H t,v [G] £ D £_1 '°°, and the statement is proved. □ 
Let us consider the following simple result. 

Lemma 5.2. Let {T n } n C V> e,oc be a sequence of r.v. 's such that forp > 1 and R > 0, 

\\ T nh, P p n < 00, \p\ < R. 

n>0 

Then the r.v. 

T( P )=Y J TnP n , \P\<R, 
n>0 

is well posed, belongs to O^ 00 and for any multi-index a with \a\ — j < t, 

D a T{p) = Y,D a T n p\ \p\<R. 

n>0 

Moreover, if 'ip £ Dom^S) then T(p)ip £ Dom^S) and 

5(T(p)V0 = ^5(T n ^)p n , \p\<R. 

n=0 
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Proof. The first assertion immediately follows by recalling that D £,p is closed w.r.t. || • \\i )P . 
As for the second one, since Dom p (5) is closed w.r.t. || ■ \\s jP , it is enough to prove that 
En=o \\ T ni>\\s, P P n < oo for \p\ < R. By using §M) one has \\T n tp\\ 5iP < &i >P ||2 1 n ||i i p||^||i i? for 
suitable p,q > p, so that 

oo oo 

W T n^hpP n < C p U\\ 5 , q W T n\kpP n < OO, \ P \ < R. 
n=0 n=0 

□ 

We can now state the key result of this section. 

Lemma 5.3. Assume that (H) and (M) hold. For k > 0, set 

A^(p) = 1 and for k > I, Af{p) = H^A^p)} + d.Af^p). (42) 

where % t,v [-\ is given in (JI) and 

r T o f T 

~<t,v _ / f/„,(,»UDl r / f/„,t,v\ j d2 



G T = ] t f( v u) dB u--r=£ J t m v )dBl (43) 

Then for every k G N, A 1 jf(p) is well posed, belongs to for all £ G N and can be represented 
as 

oo 

Af(p) = E A ^> n ' 

n=0 

with {A£ n } n given by 

• k = 0, n > : A 'q = 1 and A '^ = /or a// n > 1, 

• t > 1, n > : A,„ = _[^_«,(_^) + 

/ V t ' v A t,v m 

^ )\ + ( n+1 ^-l,n + l 



l:2l+l<n 



in which 



V* = £ f{v?)dB\ and V*> v = £ f(i$ v )d& a . 



Proof. We prove that for every k > 0, there exists a sequence {A fc '^} n independent of p such 
that 

• {AjjJJn C D A °° for every leff; 

oo 

^ (45) 



Af(p) = E A M^ iH<i- 



n=0 
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If (145]) holds, we can apply Lemma [5 .2\ so that for |p| < 1 we get A^(p) £ D f, °° for all £. And 
as a consequence of the power series expansion, A^(p) is continuously different iable w.r.t. p. 
We prove now (T4"5]) by induction. 

If /c = 0, fj4"5]) trivially holds. So, we assume that f l4"5]) is true for k — 1, A; > 1, and we 
prove that fT4"5]) holds for /c as well. For the sake of simplicity, set U t,v = f(v t s ' v )dBl and 

V*' v = jf f(vl' v )dB 2 s , so that Gf = U l ' v - p(l - p 2 )" 1 / 2 ^'. Notice that under (H) one has 
U f ' v , V f ' v £ for all i. From flU]) we have 

n [G.A^ip)}- T _ t s 2 [ f{i;tn ) 

. {1 - p) H jw) )- p{1 - p) H f{ V *>°) ) 



T-t 



We prove now that one actually has U*> v Aj^p) / / (v^) , V** tig_ x (p) / f (v*' v ) £ Dom 00 {8 2 ). So, 
we set -0 such that "0i( s ) = and ijj^s) = V/C^s")- We stress that (H) and (M) both give 
■ip £ Dom 00 (5). Since (|4"5]1 holds for fc — 1, we can write 

TTt,V\t,V / \ OO 

u ^fc-lW rrM a*. 1 ' ; / \ n I I 

/(^ ) ^ 

where, for every £ £ N, {A^ ln }„ C IT/' 00 and for any p one has Yl™=o ll^-l-i JkpP n < 00 
if |p| < 1. So, by Lemma I5T2] it is sufficient to prove that Y^=o ll^^-^k-i n \\i,pP n < 00 f° r 
|p| < 1. From pi}, there exists p,q> p such that H^'K-i.Jkp < ai iP ||f/'' 1 '||i i p HA^Ji^, so 
that 

oo oo 

E n^ A l-i,niii,pp n < E IK- Jk^" < °°- 

n=0 n=0 

This gives U t ' v A t f l^_ 1 (p)/f(v t ' v ) £ Dom oc {82) and moreover, 



5> 



Similarly, we have that V t,v A^ 1 (p) / / (v t,v ) £ Dom OQ {82) and 



So, we can write 

1 °° r /U t ' v \ t,V \ /V t ' v \ t ' V 

^ L v /(«.* ) y v /(«.' ) 

Now, since (1 - p 2 )~ 1/2 = ^ n > p 2n and p(l - p)" 1 = ^ n > p 2n+1 , straightforward 
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computations give 



n=0 l>0:2£<n 



2C 



(£\y¥ <\ f(v.) 



+ 



E * 



e>0:2£+l<n 



V ' V Ak-l, n -2l-l 



Moreover, one has 



with 



n=0 



E( rl + 1 )ll A ^ + lll^ n<0 °' IPl < 1 



n=0 



by the square-test. Therefore, by resuming we get Ajf(p) = Y^=o^knP n with 



A 



fc.n 



1 r V (21)! / ^A fc _ 1;n _ 2 A 
T-t L ^ (£!) 2 4^ 2 V fiv^) J 



£>0:2£<n 



E * 

£>0 : 2^+l<n 



/(«; 



)] + (^+i)Ai'\ n+1 , 



which gives formula (PH|) . Now, in order to conclude (and to justify the above change of order 
of the series indexes), it remains to show that 

oo 
n=0 

for all i and p. Notice that, by re- adjusting the indexes of the series, it is sufficient to prove 
that if |p| < 1, 



T\M ut,VA f: 1,n ) p n <oo and yh 2 ( 



yt,v \t,v 



n=0 



n=0 



p n < OO 



(recall that £^° =0 (n + ^IK- 

.i jn+ i\\e,pP n < oo 5 \p\ < !)• Consider the first series and, for the 
sake of simplicity, let us set g as a 2-dimensional process whose first component is null and 
its second one is given by U t,v / f (v t,v ) . Thanks to (H) and (M) one gets g e L 1 ' 00 for all i. 
Moreover, 

'U t ' v K t ' v _ \ f T 

1 = A k _ 1>n 6(g) - J (D s Al v _ 1>n ,g s )ds, 



f(v hV ) 



so that 



<l|A£V l %)lk P + (DJiWg.)ds 



< C l J\\A\\ n \\ ltP \\5(g)\\ l) q+ \\A\ V _ X 



< a 



C.p 



l|Afe!. ln ||^ + i,p, 
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for suitable p,q > p and a positive constant C^ p coming from the use of inequalities (138j) 
and (jlOjl . By recalling that, by induction, one has ^^L ll^-l-i nll^+i,p < °°' one S e ^ s 

EZo\M^j§^)\kpP n < °°- Similarly, one proves that £~ J5 2 (^^)||,, p p" < oo, 
and the statement is now completely proved. □ 

We are now ready to prove our main result. 

Theorem 5.4. Assume that (H) and (M) hold. Letip denote a Borel function with exponential 
growth, that is \4>(x)\ < C(l + e Xx ), for some C, A > 0. Then for every k > 1 one has 

d'M^n) = E(<K^) Af (p)) (46) 

where A t P v (p) is given by ( f]p| ). In particular, one has 

^E(0(^)) |p ^ o =E(0(^)Ai';o) 
where Ajj,' q i/ie starting point of the sequence {A^}„ defined through 

Proof. We first assume that G C^°, the general case to be treated by regularizing arguments. 
For k = 0, (|46i) trivially holds. By induction, we assume that (1461) holds for k — 1, k > 1, and 
we prove that (146]) holds for k as well. In fact, by recalling that d p ^ ,v = G p ,v , one has 

^E(0(^))=9 p E(0(^)A' fc ! 1 (p)) 

= E(^ w )a^Aii 1 (p))+E(^ w )a p Aii 1 (p)) 

= E ' v A£M + 0, ( P )) ) 

= E(0(^)Af(p)) 

and the above interchange between the operators E and d p is due to the fact that from Lemma 
15.31 one has A^(p) G LP for all k and p. 

Now, in order to deal with a function which is simply Borel measurable and with exponential 
growth, we use arguments similar to the ones developed in Fournie et al. [13] (see the proof of 
Proposition 3.2, at page 400) by approximating <j) with a sequence {<f) n }n C C^°. We first notice 
that <fi(£l' x ' v ) = h(Sl ,x,v ) with h having polynomial growth, so that sup^i <B sup t<s<T |0(£s X,:E )l e 
L p for every p. So, it is enough to prove that sup| p | <i? |A^'(p)| G LP for every i? < 1. And in 
fact, since 

oo 

3up|Af(p)|<^|A*fJ^ 

|Pl<R n =0 

one gets 

oo oo 

II sup |Af {p)\\\ p < Kn\\p Rn ^ E ii A £Jk^ n 

1^1^ n=0 n=0 

which is finite, as already seen in the proof of Lemma loT3"l see fj43|) . Finally, since Aj^(O) = A^,'q, 
last statement holds as well. □ 

As a particular case, useful for practical purposes, we can state the following 
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Proposition 5.5. Assume that (H) and (M) hold. Let </> denote a Borel measurable function 
with exponential growth. Then one has 

9M^t' V ))^ = E(#g*>ft) and dlE^'l)^ = E(0(g"Xo) 

where 

as; =^^((^) 2 - jf l/(^)r 2 ^) - ^ Y t,v Z t,v + 1. 

m which 

U*> v = f f{vT)dBl V* = f f{vT)dBl = f -rXsrdBl 

Jt Jt Jt j{Vs > 

Proof. We apply Theorem 15 .41 by using (|44|) . we only have to compute A{ v and v , for which 
we need also A*'". Since U t,v is independent of B 2 and Aq V = 1, we have 

Af = — !—£/*'" / T -^r- d£ s 2 = — !— t/^Z^ 

As for A*'i, by recalling that Afti = 0, (gU) gives 

A » = "fb* (7^) ) = "fb 1 7W) dB - ~ [ SS *) 

= _ v^ v Z^ v + 1 

T-t 

because D 2 s V l ' v = f(v$ v ). Finally, since D 2 S Z^ V = one has 

A t,« 1 ( ({J 1 '") 2 Z l > v \ 1 , t 

/ f T ^_ 2 r '!%z^ ds \j_ v t,v Z t,v + 1 

(T-t) 2 \ J t /(„*•«) s J t f( v *>") J T-t 

=S(^ > -jT7^ cto )-r^ W + L 

□ 

Remark 5.6. TTie formulas in Proposition 15. 51 follow from the Malliavin integration by parts 
formula in Lemma \5. 1\ which in turn is stated by means of the noise given by B 2 . This gives 
rise to simple formulas, mainly because B 2 does not appear in the dynamics for the volatility 
process. It is worth to be said that other integration by parts formulas can be provided, also 
by considering both noises B 1 and B 2 , and results similar to Proposition 15.51 can be stated. 
For example, we have found several weights that involve all the noises and lead again to an 
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expression for the derivatives of order one and two, evaluated in p = 0, written in terms of 
some suitable weights. But such resulting weights are really complicated to write down and then 
unfeasible to be used in practice. However, as we will see in next Section 0, for numerical 
purposes the weights in Proposition 15.51 work very efficiently, and this is a further reason to 
drop complications. 

Remark 5.7. The formulas for the weights continue to hold if jumps are suitably inserted in 
the model. For example, Aids, Leon and Vives ^ consider a jump- diffusion model for the log- 
returns by inserting in the dynamics |3|) an independent compount Poisson process (actually, 
in that paper the model for the volatility is more general - it is not a diffusion process, neither 
Markovian, but it is measurable w.r.t. the noise B l ) and they prove that a Hull and White type 
representation for the call price again holds (see Theorem 4-1 therein). In our framework, one 
would have 

is = x+ I (r-m\- l -f\v u )yu+ I f(v u ) (pdB l u + y/l - p* dB 2 u ) + ^ A n , (47) 

^* n=l 

where v satisfies (CP, N stands for a Poisson process with intensity A > 0, the sequence of 
the jumps {A n } n is modeled by i.i.d. r.v's and the random sources B, N and {A n } n are 
independent. Since we consider the evolution of the log-returns under the risk neutral measure, 
in Q?U one has m = J (e z — l)p(dz) , p denoting the common law of the A n 's. Now, if we assume 
that e Al G LP for every p (and this holds in cases of interest, e.g. when e Al follows a log-normal 
or an exponential law), then it is easy to see that Proposition 15.51 continues to hold with £ 
fulfilling as a consequence of the fact that the noise from the jumps is independent of the 
Brownian noise. This is done e.g. in Bally, Bavouzet-Morel and Messaud JE/, where a Malliavin 
calculus in the direction of the jumps is also developed, possibly giving other integration by parts 
formulas. 

Remark 5.8. In practice, it is well known that the Malliavin representation for sensitivities 
leads to a problem of high variance. A standard variance reduction technique uses localizing 
functions, as follows. Let h denote a pay off- function, assumed to be Borel measurable and with 
polynomial growth. By Proposition \5.5[ for a suitable regular and with bounded derivatives 
function $ ; we can write 

dpEQiiSp"")) = - ^)(S^ x ' v )A{ v (p)) +E(d x ^S^ v )d p S t /' v ) 

and similarly, 

d 2 p E(h(S^ v )) = E((h - $)(S^)Af (p)) + E(dmS^ x ' v )(d p S^ x ' v ) 2 + d^iS^dlStf'*). 

So, the idea is to take $ such that h — $ cannot assume values too different each other. This 
procedure permits to decrease the contribution given by high fluctuations of the weights A*' w (p) 
and A^p) and, as a consequence, it improves the numerical results and reduces the variance. 
The literature on these arguments already contain many different possibilities for the localizing 
function We need here $ G C 2 that allows to work efficiently for the second order repre- 
sentation. So, in our numerical analysis we consider a slightly modification (more regularity) 
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of the localizing function for call options developed in Bavouzet- Morel and Messaud JPJ/. In 
particular we take 




— oo 



w 



B$(z)dz dw with B$(z) 











(z - Kf 3 



45 3 A8 



ifz<K-5 

ifze[K-5,K + 5] (48) 
if z > K + 5 



6 Numerical examples 

We study here the numerical behavior of the series coefficient representations: we use the results 
of sections H] and [5] to numerically compute by means of a Monte Carlo method the call price 
through its series expansions up to order 1 and 2 and we compare the results with the ones in 
Antonelli and Scarlatti [5] (where it is clear that only the terms up to order 2 are numerically 
relevant). 

The tests in [5] are related to the most popular stochastic volatility models: 
• the Hull and White model [17], in which 



These are models for which the validity of the requirements (H), (A) and (M) deserves some 
comments (we recall that (AS) always holds, as proved in [5], Appendix A). First of all, the 
fact that (H) is not fulfilled for the Heston model is not really a problem (as already observed 
in [5], Section 5.1). Concerning (A) and (M), they are not generally verified by the models we 
are concerned in. For example, (A) always fails. This is not really a problem in the Hull and 
White model, since a transformation can be done in order to get f(v) = const > 0, but it does 
become an unpleasant fact for the other models. However, in order that (A) and therefore 
(M) hold, one could think to perturb a little bit the model. Since problems come from a 
function of the type g(v) = v a , a = 1, 1/2, one could replace g with g e (v) = \/v 2a + e. Such a 
perturbation procedure has been exploited in our numerical simulations concerning the Stein 
and Stein model and the Heston model. In any case, we have set e = 10~ 5 . 

Let us now enter in the details of the simulation procedures. 

We use an Euler approximation scheme to simulate the dynamics of the diffusions over [0,T] 
(in all tests, we have set t = 0). For the simulation of the volatility process in the Heston model 
(following a CIR process), we have done a further approximation: in the approximating scheme, 
the diffusion coefficient n(y) = y/v of the volatility process is replaced by rj y (v) = \/\v\ + 7. In 
practice, we have set 7 = 10~ 5 . It is worth to be said that the use of more accurate schemes for 



fJ>(v) = fiv, T)(v) = cv, f(v)=v; 



• the Stein and Stein model [12], in which 



fi(v) = b(a-v), r](v) = c, f(v) = v; 



• the Heston model [TB], in which 



H(v) = b(a-V), 7](v)=Cy/v i f(v) = y/v. 
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the simulation of the CIR process could be used, thus avoiding this further complication. For 
example, Alfonsi [T] has developed a second order approximation scheme for the Heston model, 
but it deserves to be modified in order to simulate the functionals we need in our formulas. 
However, as it follows from the numerical tests, our simulation scheme behaves efficiently. 
So, in our experiments we use the above Euler approximation scheme to simulate the dynamics 
of the diffusion processes; moreover, a Monte Carlo algorithm is developed to compute the 
expectations giving the coefficients of the power series. We divided the time interval [0,T] in 
N = 500 parts and we have performed 10 4 simulations, and this means that we have used a 
standard number of parameters for the Euler scheme and for the Monte Carlo replications. 
The numerical results refer to the percentage errors from the different expansions up to order 
1 and 2. For the sake of clearness, we precise that the "percentage error" is given by 

1^1 x 100. 
V 

where p denotes the true option price (or a benchmark value for it) and p stands for the 
approximation for p. 

The tables give the percentage errors in our numerical experiments from the various methods 
and are labeled as follows: 

• 'AS' displays the results in Antonelli and Scarlatti [5]; 

• 'ExpA-1' and 'ExpA-2' refer to the Taylor series expansion up to the first and the second 
order respectively, obtained from the results developed in Section HJ 

• 'ExpM-1' and 'ExpM-2' relate to the Taylor series expansion up to the first and to the 
second order respectively, obtained in Section [5] and refined with the localizing function 
(T48]h in which 5 = lO -2 ^, K denoting the (varying) strike price. 

We stress that the results labeled 'AS' come from the original paper [5] (see Section 5 therein) 
and we recall that they are computed through the first order expansion by means of a further 
approximation for the coefficients go and g\ (see d§J) and (ITUl) respectively): go and g\ are 
changed with g and g x obtained by replacing (M t,v )[ t)T } with F,((M t,v )[ t> T\). It is clear that 
~g Q and g 1 give a ready first-order closed-formula approximation for the call price, which is 
very accurate. So, the method 'ExpA-1' gives the first order expansion from the Antonelli and 
Scarlatti approach deprived of any further approximation, and we will see that is works better 
in several cases. 

The tables are collected in Section [71 They refer to varying values of p and of the strike price 
K. As for the other parameters, the choice is done according to the model and can be detailed 
as follows. 

• Hull and White model - Table [TJ 

As in |5], we have set: t = 0, T = 0.5, S = 100, v = 0.2, r = 0.0953, c = 0.1 and 
/i = 0.2. Table [T] shows different percentage errors as p and K vary. The percentage error 
is computed with respect to a benchmark value, obtained by evaluating the price of the 
European call option with a Monte Carlo method when the number of iteration is 10 6 
and the time interval is divided in 10 3 parts. 
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• Stein and Stein model - Table [2j 

Again, the parameters are set as in [5], that is: t — 0, T — 0.5, So = 100, vq = 0.2, r = 
0.0953, a = 0.2, b = 4, c = 0.1. The percentage errors are eveluated w.r.t. the exact 
price. 

• Heston model - Table [3] and HI 

The square-root model for the volatility process ensures that v a.s. never reaches if 
v o > and the condition 2ab > c 2 (often called Novikov condition) is fulfilled, and in 
such a case one has also nice Malliavin derivability properties for v (see Alos and Ewald 
[3]). The results in [5] are given when the Novikov condition does not hold. So, we present 
for this model two sets of results: 

— Table [3] refers to parameters chosen in Alos [2]: t = 0, a = 0.04, b = 8, c = 0.1, 
r = 0.0953, v = 0.0225, S = 100, so that condition lab > c 2 is fulfilled; 

— Table H] refers to parameters chosen in Antonelli and Scarlatti [5j: for short term 
options (T = 0.1 or T = 0.2) then p = -0.76, a = 0.025, b = 1.62 c = 0.44; for 
long-term ones (T = 0.4, T = 0.5 or T = 0.8) then p = -0.64, a = 0.035, b = 1.15, 
c = 0.39. Here, all cases give 2ab < c 2 . 

The errors are computed w.r.t. the exact price, evaluated by using the pricer at http:// 
kluge . in-chemnitz . de/tools/pricer/. 

It appears evident that Expansion A gives highly efficient results. In fact, the 'ExpA-1' and 
'ExpA-2' percentage errors are always of order strictly smaller than 1%. In addition, since 
'ExpA-1' coincides with the Antonelli and Scarlatti expansion up to the first order, it is also 
clear that often the approximations ~g Q and g 1; giving a closed-form first order approximated 
formula, bring to loose accuracy. Moreover, the use of the second order coefficient generally 
improves the results (as shown in the 'ExpA-2' errors), sometimes cutting the error of some 
order, and this is not really a problem from the numerical implementation point of view: the 
ExpA algorithm takes great advantage from the fact that the structures of the coefficients 
of order 1 and 2 are very similar each other, while the implementation of the second order 
coefficient from the Antonelli and Scarlatti expansion is actually heavy. 

Concerning Expansion M, i.e. the one from the Malliavin approach, in the Hull and White 
model and in the Stein and Stein one the errors do not exceed 2%, so perfectly in the desired 
range. Such a performance also holds for the Heston model when the Novikov condition is 
fulfilled (Table [3]). But if such a condition does not hold (Table H]), the errors generally increase, 
even if they remain below the standard 5% threshold but with two exceptions (ExpM-2: p = 
-0.64, K = 110 and ExpM-2: p = -0.76, K = 110), in which they are of order 6%. This 
could be explained by recalling that when the Novikov condition is not satisfied, the Malliavin 
differentiability of the volatility process is not guaranteed (see Alos and Ewald [3]). 
Nevertheless, our opinion is that the largely quite accurate behavior of the Malliavin method 
is highly satisfying, mainly for two reasons. Firstly, as already observed, it is ready to be 
implemented for more general European options (even with a non regular payoff-function) and 
moreover, following Remark 15.7} it is worth to be used also if the log-returns are driven by the 
more general jump- diffusion model in (14T|) . 
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7 Tables 





K = 90 


K = 95 


K = 100 


K = 105 


K = 110 


p = -0.25 


AS 


0.386235 


0.440003 


0.435233 


0.339402 


0.127030 


ExpA-1 
ExpA-2 


0.263909 
0.264210 


0.309145 
0.309780 


0.353728 
0.354730 


0.395341 
0.396792 


0.438145 
0.440186 


ExpM-1 
ExpM-2 


0.199430 
0.284079 


0.425136 
0.496970 


0.469145 
0.319980 


0.725978 
0.267020 


0.910577 
0.514088 


p = -0.5 


AS 


0.169833 


0.141649 


0.002633 


0.335407 


0.940766 


ExpA-1 
ExpA-2 


0.041942 
0.040742 


0.080049 
0.077808 


0.117349 
0.114099 


0.152940 
0.148380 


0.196648 
0.189953 


ExpM-1 
ExpM-2 


0.007453 
0.415920 


0.124952 
0.529802 


0.566645 
0.795568 


0.239341 
0.518876 


0.315632 
1.636366 



Table 1: Hull and White model - percentage errors w.r.t. a benchmark value. Parameters: t 
T = 0.5, S = 100, v = 0.2, r = 0.0953, c = 0.1, n = 0.2. 





X = 90 


K = 95 


K = 100 


K = 105 


K = 110 


p = -0.25 


AS 


0.003941 


0.123226 


0.326447 


0.583335 


0.797076 


ExpA-1 
ExpA-2 


0.007070 
0.005262 


0.003160 
0.000415 


0.011961 
0.014921 


0.032119 
0.035616 


0.039606 
0.046035 


ExpM-1 
ExpM-2 


0.258058 
0.163152 


0.240485 
0.027884 


0.231560 
0.069982 


0.167886 
0.563581 


0.149144 
1.689356 


p = -0.5 


AS 


0.047087 


0.160828 


0.332411 


0.571909 


0.854295 


ExpA-1 
ExpA-2 


0.051959 
0.044712 


0.057347 
0.045937 


0.024588 
0.011695 


0.012997 
0.028809 


0.035652 
0.007137 


ExpM-1 
ExpM-2 


0.336514 
0.295035 


0.383521 
0.307504 


0.449586 
0.506819 


0.293623 
0.527940 


0.452388 
0.521331 


p = -0.75 


AS 


0.144578 


0.239578 


0.340785 


0.522829 


0.950355 


ExpA-1 
ExpA-2 


0.128902 
0.112629 


0.157840 
0.130979 


0.133673 
0.101239 


0.106995 
0.065219 


0.285850 
0.212845 


ExpM-1 
ExpM-2 


0.390693 
1.650679 


0.164782 
0.136951 


1.133769 
1.524349 


1.812543 
1.675562 


1.634112 
1.492334 



Table 2: Stein and Stein model - percentage errors w.r.t. the exact price. Parameters: t = 0, T 
S = 100, t; = 0.2, r = 0.0953, a = 0.2, b = 4, c = 0.1. 
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K = 90 


K = 95 


K = 100 


K = 105 


K = 110 


p = -0.5 


T = 0.5 


ExpA-1 
Exp A- 2 


0.006135 
0.005404 


0.007329 
0.005785 


0.003350 
0.001136 


0.002555 
0.005671 


0.002873 
0.007891 


ExpM-1 
ExpM-2 


0.128272 
0.663177 


1 22593 
0.431397 


068726 
1.190316 


1 35990 
0.512402 


060668 
0.659557 


T = 0.8 


ExpA-1 
ExpA-2 


0.001454 
0.000532 


0.001302 
0.002577 


0.008102 
0.009491 


0.018863 
0.020320 


0.031219 
0.033123 


ExpM-1 
ExpM-2 


0.434218 
0.681888 


0.500039 
0.145416 


0.605644 
0.853767 


0.373685 
0.912884 


0.142499 
0.473907 


T = 1 


ExpA-1 
ExpA-2 


0.021982 
0.020960 


0.031822 
0.029948 


0.042954 
0.040171 


0.055580 
0.051814 


0.071431 
0.066535 


ExpM-1 
ExpM-2 


0.388991 
0.351164 


0.553972 
1.212113 


0.872142 
0.845210 


0.734878 
0.171887 


0.642501 
1.006377 



Table 3: Heston model with parameters fulfilling the Novikov condition - percentage errors w.r.t. the 
exact price. Parameters: t = 0, S = 100, v = 0.02225, r = 0.0953, a = 0.04, b = 8, c = 0.1. 





K = 90 


K = 95 


K = 100 


K = 105 


K = 110 


p = -0.64 


T = 0.4 


AS 


0.004745 


0.052284 


0.096813 


0.149266 


0.231650 


ExpA-1 
ExpA-2 


0.145056 
0.076859 


0.084289 
0.014733 


0.033238 
0.041668 


0.052816 
0.040969 


0.213478 
0.078879 


ExpM-1 
ExpM-2 


0.049668 
1.682665 


0.169833 
2.722891 


0.040386 
1.967607 


0.330853 
1.843184 


0.070622 
3.016751 


T = 0.5 


AS 


0.004457 


0.058372 


0.109874 


0.170933 


0.264214 


ExpA-1 
ExpA-2 


0.217793 
0.121799 


0.163739 
0.065143 


0.124153 
0.018367 


0.159125 
0.031568 


0.339199 
0.166553 


ExpM-1 
ExpM-2 


0.295190 
0.956495 


0.250897 
0.738435 


0.348796 
3.272917 


0.545558 
1.394025 


0.301669 
6.177427 


T = 0.8 


AS 


0.023857 


0.046118 


0.115713 


0.201686 


0.326068 


ExpA-1 
ExpA-2 


0.040934 
0.134341 


0.094668 
0.274656 


0.230403 
0.420120 


0.316294 
0.531201 


0.292525 
0.557421 


ExpM-1 
ExpM-2 


0.269246 
2.030071 


0.629545 
0.960147 


1.258894 
2.772005 


1.664970 
1.073997 


2.582807 
3.222600 


p = -0.76 


T = 0.1 


AS 


0.044528 


0.008180 


0.051279 


0.057639 


0.051020 


ExpA-1 
ExpA-2 


0.075473 
0.061353 


0.031628 
0.016046 


0.034223 
0.051183 


0.036352 
0.002745 


0.479536 
0.399934 


ExpM-1 
ExpM-2 


0.715144 
1.690227 


0.086117 
4.377866 


0.484912 
0.592030 


0.177757 
4.588024 


0.709826 
5.912330 


T = 0.2 


AS 


0.043600 


0.029069 


0.090594 


0.121116 


0.082985 


ExpA-1 
ExpA-2 


0.138635 
0.094970 


0.057408 
0.012563 


0.030773 
0.079821 


0.014847 
0.059268 


0.395983 
0.259167 


ExpM-1 
ExpM-2 


0.390806 
1.712282 


0.121161 
0.612300 


0.747368 
1.215959 


1.999056 
1.571722 


1.871387 
2.452193 



Table 4: Heston model with parameters not fulfilling the Novikov condition - percentage errors w.r.t. 
the exact price. Parameters: t = 0, S = 100, v = 0.02225, r = 0.0953; for p = -0.64: a = 0.025, 
b = 1.62, c = 0.44; for p = -0.76: a = 0.035, b = 1.15, c = 0.39. 
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